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Abstract 

We show that the Fade Approximant (PA) approach for resummation of perturbative se¬ 
ries in QCD provides a systematic method for approximating the flow of momentum in 
Feynman diagrams. In the large-/3o limit, diagonal PA’s generalize the Brodsky-Lepage- 
Mackenzie (BLM) scale-setting method to higher orders in a renormalization scale- and 
scheme-invariant manner, using multiple scales that represent Neubert’s concept of the dis¬ 
tribution of momentum flow through a virtual gluon. If the distribution is non-negative, 
the PA’s have only real roots, and approximate the distribution function by a sum of 
5-functions, whose locations and weights are identical to the optimal choice provided by 
the Gaussian quadrature method for numerical integration. We show how the first few 
coefficients in a perturbative series can set rigorous bounds on the all-order momentum dis¬ 
tribution function, if it is positive. We illustrate the method with the vacuum polarization 
function and the Bjorken sum rule computed in the large-/3o limit. 
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1 Introduction 


Fade Approximants (PA’s) are known to be useful in many physics applications, 
including quantum field theory and statistical physics 0] . These applications of the 
PA method have recently been extended to QCD |^, ^], where the method has been 
shown to be effective both in predicting unknown higher-order coefficients and in 
summing the perturbative series. In these applications of PA’s to QCD, the generic 
starting point is a perturbative series for some physical observable A that has been 
calculated exactly to some hnite order 

A ~ An = x(^Cq CiX C 2 + ■ ■ ■ CnX^) (1) 


where the the expansion parameter x is related to the renormalized coupling constant 
by X = x(/i^) = /{4:7i), where fi is the renormalization scale in some scheme 

such as MS . The corresponding PA’s are ratios of polynomials 


x[N/M] 


9 /V 

Po -h pix + P2X ... pnx 
^ _ 

1 qix q2x‘^ + ... qMx’^ 


( 2 ) 


with N + M = n, chosen such that they reproduce the known coefficients Co through 
Cn when expanded back in a Taylor series. It is clear that a PA (|^) includes some 
higher-order effects: the hope is that it resums physically relevant higher-order 
contributions, so that x[N/M] will be a better approximation to the observable A 
than the original truncated power series An- 

In the absence of exact perturbative calculations for QCD observables beyond 
three loops, the evidence for the effectiveness of the PA method in QCD applications 
is mostly indirect: 


• As already mentioned, PA’s are successful in resumming the series as well as 
in predicting higher-order coefficients in other models 1,0: see also 0 and a 
recent review in [^. Among these, an important example @] is provided by the 
limit of QCD where the number of light flavours Nf becomes very large, al¬ 
though there are important differences in the physics described by this theory, 
due to its lack of asymptotic freedom. The large-order renormalon behavior of 
the perturbative coefficients is thought to resemble that of QCD, and there is 
a theorem that PA predictions for higher-order perturbative coefficients must 
converge if the perturbative series is dominated by a renormalon. 

• Comparisons of the PA method p| with other methods that seek to optimize 
the perturbative result through a proper choice of scale and scheme, such as the 
method of Effective Charges 0 , the Principle of Minimal Sensitivity, 0] and 
the BLM scale-setting method 0 show good numerical agreement for the 
Bjorken sum rule and in certain cases also exhibit close algebraic relations 0. 

• PA’s were found |Q] to reduce the undesirable renormalization scale- and 
scheme-dependence of physical observables, as compared to the partial sums 
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on which they are based. Some understanding of this result was provided in |P] , 
where it was proven that diagonal x[N — 1/A^] PA’s become exactly scale in¬ 
variant when the (3 function is approximated by its leading term. This strongly 
suggests that PA’s resum correctly higher-order contributions associated with 
the running of the coupling constant: see also 0 for recent intriguing work 
along a somewhat related approach. 


Despite these pieces of evidence for the relevance of PA’s for QCD, there has so 
far, to our knowledge, been no direct diagrammatic interpretation of the summation 
by PA’s. The magnitude of this problem can perhaps be seen from the fact that a 
Pade function has simple poles in the coupling-constant plane, and therefore can¬ 
not reproduce the expected factorial growth of the perturbative coefficients. This 
factorial growth is apparently an essential feature of higher-order contributions, ap¬ 
pearing because of the flow of extremely high or extremely low momentum through 
virtual gluon lines, ultraviolet (UV) and infrared (IR) renormalons, respectively, and 
also because of the multiplicity of higher-order diagrams. 

It is clear from these considerations that PA’s cannot account for the full set 
of higher-order graphs. On the other hand, as reviewed above, there are strong 
indications for the relevance of the summation of higher orders of perturbative QCD 
by PA’s. How can we interpret this summation in terms of Feynman diagrams? 

Studies of higher-order perturbative QCD diagrams are often made by hrst de¬ 
composing them in a skeleton expansion, in which each term contains chains of vac¬ 
uum polarization bubbles inserted in virtual-gluon propagators. These have been 
studied in the BLM approach, which seeks the optimal scale for evaluating each 
term in the skeleton expansion. The last step, the sum over skeleton graphs, is then 
similar to summation of perturbative contributions for a corresponding theory with 
/? = 0, i.e., a conformal theory |^, |^. We shall adopt a similar procedure here. 

In this paper, we consider a subset of graphs corresponding to a single virtual- 
gluon exchange. We adopt the concept of the momentum distribution function intro¬ 
duced by Neubert (see also [|^]): in the large-/do limit, the all-order summation 
of diagrams is reduced to a single integral over all scales of the running coupling 
constant, with a weight function describing the distribution of momentum flowing 
through the gluon line [|. This is natural in QED calculations where the standard 
running coupling a{—k‘^) sums all vacuum polarization corrections to the photon 
propagator. In QCD, the same feature is incorporated into the scheme 

dehned from the potential for the color-singlet scattering of two heavy colored test 
charges. Since the coupling is singular in the large-/?o limit, the integration over 
the gluon momentum yields in general renormalon singularities. We hnd that PA’s 
make a systematic approximation to the momentum distribution function, thereby 
extending the leading-order BLM prescription. 

In many physical cases, the distribution function has been found empirically to be 


Mhe relation of the optimal renormalization scale and the BLM prescription to a weighted 
momentum flow integral is also discussed in refs. [11^ Iff 
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non-negative If the distribution function is indeed non-negative, the resummed 
amplitude dehnes a so-called Hamburger function Under this assumption, we 
obtain the following results: 


• One may interpret the x[N — 1/A^] PA as a discrete approximation to the 
integral. This is because, for Hamburger functions, the poles of the diagonal 
PA’s are all real, corresponding to meaningful physical scales, and their weights 
are all positive. 


Moreover, there is a formal connection between PA’s and the Gaussian quadra¬ 
ture method for numerical integration ||^. The basis for this connection is the 


observation that using PA’s amounts to approximating the distribution func¬ 
tion by a sum of 6 functions. The locations of these S functions are determined 
by the poles of the PA, and their weights are determined by the residues. 


• Furthermore, one can use a known PA-based method to bound rigorously the 
all-order distribution function, by using only the hrst few coefficients of the 
perturbative series. 


It is natural to ask in addition whether higher-order PA’s converge to the re- 
summed result, i.e., to the value of the integral over the exact, continuous distri¬ 
bution function. This question is well posed only when the integral itself is well 
dehned. However, it is well known that this is not the case in the presence of 
the IR renormalons expected in QCD, and indeed we hnd that in general PA’s do 
not converge []. Only if there are no IR renormalons, meaning that the momen¬ 
tum distribution function completely vanishes for “IR scales” in addition to being 
non-negative, in which case we would have a Stieltjes function rather than just a 
Hamburger function, would the integral be well dehned on the positive real axis in 
the as plane. In this case, higher-order PA’s do converge to the correct value. 

The reader may hnd it useful if we relate our approach to that of [§] . This is based 
on choosing the renormalization scale of the coupling constant so that the next-to- 
leading (NLO) coefficient in the leading /3o series vanishes. The physical meaning of 
this scale setting is that certain higher-order vacuum-polarization ehects are taken 
into account. The analysis of |jT^ is a generalization of the BLM approach. The 
BLM choice of scale amounts to approximating the continuous distribution function 
p{s) by a single S function IjT^: p(s) —^ ^0(5(8 — Sblm)- As shown in @], 


m 


the large-/3o limit, the leading-order BLM procedure is also exactly equivalent to a 
x[0/l] PA, i.e., the BLM procedure is equivalent to a single geometrical series in the 
renormalized coupling constant at an arbitrary scale. Here we show that higher-order 
diagonal a;[A^ — 1/A^] PA’s in the large /do limit can also be described naturally in the 
language of momentum distribution functions: they correspond to approximating 


Hn some cases the UV cut-off must be chosen appropriately in order to achieve this. 

§See and section 3.1 for the exact definition. 

^Although the PA predictions for the next individual terms in the perturbation series become 
progressively more accurate. 
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the distribution function by N ^-functions, a multi-scale extension of the BLM idea. 
Several interesting extensions of the BLM method have been suggested in the 
past [T^, [T^, [^. The motivation was to include higher-order effects, both within the 
large-/?o approximation and outside it. These suggestions were based on improving 
the leading-order BLM scale [|T^ |T^, on introducing multiple scales to account for 
non-leading Pq terms [^, improving the single BLM scale at higher order |[T6[| , 
and on setting the scheme by tuning the coefficients of the (3 function [^. None 
of these extensions, however, introduces multiple scales already within the large-/3o 
approximation, and in this sense these methods have no natural description in the 
physically attractive language of momentum distribution functions, as is done by 
diagonal PA’s. 

The outline of the paper is as follows: in the next section we review the BLM 
approach and the concept of a momentum distribution function. In section 3 we 
study the relation between PA’s and approximations to the momentum distribution 
consisting of a sum of 5 functions. We use some known mathematical results con¬ 
cerning Hamburger functions to draw the above-mentioned conclusions concerning 
the use of PA’s in QCD. In section 4 we illustrate the application of our method to 
the vacuum polarization D function and the Bjorken sum rule in the large-/do limit. 
Finally, in section 5 we give our conclusions and discuss some questions raised by 
this work. 


The BLM Method and the Momentum Distri¬ 
bution Function 


2.1 Definition of the Momentum Distribution Function 

The BLM method seeks to absorb most of the growth of the higher-order pertur¬ 
bative QCD coefficients with an astute choice of renormalization scale in lower-order 
expressions. Formally, one arranges the perturbative series for a generic QCD ob¬ 
servable in a skeleton expansion, whose coefficients are given by a conformal theory, 
and then seeks an optimal scale to evaluate each term in the skeleton expansion. 
The relation of this approach to the a;[0/l] PA in the large-/do approximation has 
been discussed in p. In order to go further, it is convenient to adopt the approach 
of 0, which considers the resummation of all orders of perturbative corrections 
to physical observables that depend on a single external momentum Q^. The re¬ 
summation method described takes into account only the subset of graphs that can 
be described as an exchange of one effective virtual gluon, which is the simplest 
example of a term in the skeleton expansion. Resummation is achieved by using 
the running coupling constant as{—k‘^) at the vertices, where k is the momentum 
flowing through the virtual gluon, instead of a renormalized coupling at some 

hxed scale jn. 
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The result of the improved calculation is denoted by Ares- 

/ fOO 

d'^kg(k,Q^) as(—k^) = J dtw^{t) x^{tQ‘^) (3) 

where g{k, Q"^) is the integrand of the Feynman diagrams. Clearly, an infinite number 
of diagrams is taken into accounted by this integral. The function w^{t) in (|]) is 
interpreted as the momentum distribution function characterizing the virtuality of 
the exchanged gluon. It provides the weighting in the integral over the running 
coupling x^{tQ^) = af{tQf) /(47r). The superscript R stands for the renormalization 
scheme, and serves to remind us that both the distribution function and the running 
coupling depend on the scheme, whilst the resummed result Ares should be scheme- 
invariant. 

In practice, for physical examples. Ares has been calculated to all orders only 
in the large-/3o approximation, i.e., when the running of the coupling is controlled 
entirely by the 1-loop (3 function. Then 


x^{tQ^) 


1 -h ln(f)/9oX^(Q^) 


( 4 ) 


where /3o = y-^c — = 11 — |iV/ for QCD. All-order calculations in this ap¬ 

proximation are possible because the perturbative coefficients are proportional to 
the coefficients of the large-A^/ limit [^, ^ in a non-Abelian theory such as QCD, 
in which the higher-order corrections are only due to fermion-loop insertions in the 
gluon propagator. It has been argued [^, that, since (3q is relatively large in 
QCD with a few light flavors, this approximation should be quite good, and we 
adopt it here. 

We note that the integration in (^) includes both high momenta: f —>■ cxo in 
the deep UV region and low momenta: f —> 0 in the deep IR region. Whilst the 
UV integration region is well defined, since renormalizability guarantees that w^(t) 
decreases fast enough, the IR integration region is ill defined, due to the Landau pole 
present in (^). This is how IR renormalons appear in this formulation, reminding 
us that perturbation theory in QCD is not adequate to describe well long-distance 
physics with a truncated lowest-order (3 function. 


2.2 Renormalization-Scheme Independence 

In the large-/3o approximation, the dependence on the renormalization scheme can 
be removed from (H) |OT . This is because the dependence on the scheme enters only 


through the parameter C, related to the finite part of a renormalized fermion-loop 
insertion in the gluon propagator: in the MS scheme C = —5/3, and in the V- 
scheme C = 0. One observes that w^{t) depends on C, Q^, and /i^ only through the 
combination /{Q'^e^), whilst w^{t)dt does not depend on C at all. We follow [0 
in defining 


- ^ 


T = t 




( 5 ) 
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and a scheme-invariant momentnm distribntion 


w{t) : w{T)dT = w^{t)dt ( 6 ) 

In this notation, (|^) becomes 

roo 

Ares = / w{T)x{Te^Q‘^)dT (7) 

Jo 


where a;(re‘'Q^) is scheme-invariant. Mathematically, the parameter C plays here 
the same role as the scale /i, and therefore there is only one free parameter in the 
renormalization gronp. Finally, we conclnde that as long as we nse the integral 
representation of Ares, the choice of renormalization scheme does not change either 
w(t) or the hnal resnlt for Ares- We shall see shortly that this is not trne when 
one nses a hnite-order Taylor series An as an approximation for Ares, where scale 
and scheme dependence do appear. However, the nse of diagonal a:[iV — 1/iV] PA’s 
eliminates the scheme and scale dependence from the hnite-order approximation for 

i- 

For onr purposes, it is convenient to rewrite ( 0 ) again, using s = ln(r) as the 
integration variable: 

/ OO 

p(s)a:(e®+^Q^)ds (8) 

-OO 

where p{s) = tc(r = e®) e^. We then see from that 


a:(e"+^Q") = 


X 


R 




X 


R 


m 


l + {s + C- ln{fiyQ^))(3ox^{p^) 1 + (s + C)/?ox^(Q 2 ) 


(9) 


and 


or 


A — 


P(s) 


X^{jJ?‘ 


A — 

■^rp.A — 


1 -|- (s -|- C 
)0 

P(s)- 


- ln(/u^/Q^)) /Jox^(/u^ 

x^(Q^) 


-ds 


-ds 


( 10 ) 


( 11 ) 


1 + {s + C) f3ox^{Q'^) 

We can think of s as a scale parameter: the long-distance physics is described by 
p{s) for negative s, while short-distance physics is described by p{s) for positive s. 


2.3 Perturbative Expansion 

We now turn to the perturbative treatment of Ares- We write the n-th order Taylor 
expansion for Ares as 

n , 

Ares ~ = X^{p?) C^{p^) (poX^ip?)) (12) 

k=0 

where the coefficients are moments of the distribution function p(s): 

/ OO , . u 

^p(s) (s+ C - ln(/iVQ^)) ds (13) 
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The first coefficient Co is just the integral over the momentum distribution function 
p{s), i.e., its normalization, and does not depend on the scale and scheme. The 
second coefficient is related to the average momentum flowing through the 

virtual gluon, and depends on the scale and scheme through 

Cpp.^) = CpQ^) + C„\n(p.yQ^) (14) 

and 

= CpQ^) - Co (c- Iti(liW)) (15) 

where the superscript V stands for the V scheme: C = 0, and R stands for a generic 
scheme characterized by some other value of C. Higher-order coefficients are given 
by higher moments of p(s), and depend on the scale and scheme through higher 
powers of (C — ln(p^/Q^)). 

Consequently, at any finite order n there is some residual dependence of the 
partial sum on the renormalization scheme and scale, through the combination 
(C — ln(p^/Q^)). This dependence is formally of the next order in the coupling, but 
in practice it can be quite large and an inappropriate choice of the scale and scheme 
can lead to misleading results, as we show later for the example of the vacuum 
polarization D function - see Fig. 1. 


2.4 The BLM Prescription as a Narrow-Width Approxima¬ 
tion 

As was shown in and mentioned in the Introduction, the x[N — 1/A^] PA based 
on An n = 2N — 1 eliminate the scale and scheme dependence completely from 
the finite-order approximation to A^es- As was also mentioned in the Introduction, 
the BLM prescription is based on choosing a renormalization scale so that the NLO 
coefficient vanishes in the large-/3o approximation in which we work. We see from 


(|T^ that this translates into 


p{s) s + C - ln(p|i^/g2) 


ds = 0 


(16) 


and therefore 


exp C + 


JZoPis)sds 


= Q 


2p-C^{Q^)ICo 


(17) 


By construction in the BLM method, the NLO approximation A 2 to A is equal to 
the leading-order approximation Ai. We obtain from (1^) and (|l7|) 


Are. ~ Ai{BLM) = A2{BLM) = Cox^ (18) 

Note that the leading-order BLM result is scheme-invariant, whilst pblm is scheme- 
dependent. 
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Another way to obtain the BLM result, which will be important for our later 
generalization, is the following. Starting with (^, when only the hrst two coefficients 
Co and are known, one can ask for what approximation to p(s) is the integral 

representation for Ares equal to Ai{BLM). The answer is [10 


p{s) = C„i(s + C + (Cf(<3")/Co)) 


(19) 


Note that this p{s) is scheme-invariant, as it should be. The BLM method can be 
thought of as an approximation to the distribution function by a single 6 function 
located at the BLM scale. As stressed in [ITj, this is a good approximation only if 


the distribution function is narrow. Information about the width and shape of the 
distribution function are encoded in higher-order coefficients, which are higher-order 
moments of the distribution function. 

To see this in another way, we substitute (|T^ into (0), obtaining 


Ares 


Ai{BLM) = Co 




1 - {C’‘{Q^)/Co)PoxR(Q-^) 


( 20 ) 


which is another, exactly equivalent, representation of ([18|) . We stress again that 
the result is exactly scale- and scheme-invariant. We could, just as well, write it as 


•Ares ~ 


Ai(BLM) = Co 




1 - (Cf(pV/C„)/)ox«(p^)' 


( 21 ) 


In our notation, (pT]) is an a; [0/1] PA. The leading-order BLM approximation in 
the large-/do approximation is therefore equivalent to the assumption that, when 
a given perturbative series is known up to NLO: a;(Co -|- Cix), an improvement 
is achieved by continuing it as a geometrical series with the ratios of coefficients 
taken as Ci/Co, obtaining: a:(Co -|- Cix -|- C\ /Ci^x^ -|- ...). In the case of a positive 
distribution function p(s), this is probably a better approximation than using the 
NLO truncated series at an arbitrary scale and scheme. However, it is not a good 
approximation unless the distribution function is narrow. In the next section, we 
will see how a:[A^ — 1/iV] PA’s provide better approximations at higher orders, whilst 
maintaining the scale-independence property. 


Positive Momentum Distribution Functions and 
PA’s 


3.1 Positivity and Hamburger Functions 

A probabilistic interpretation of the momentum distribution function has been pro¬ 


posed in 0: clearly, an important issue for this interpretation is whether the 
distribution function p(s) is non-negative. Moreover, if the sign of p(s) changes, 
and in particular if there are large cancellations within the first moments, the BLM 
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method may give bad predictions. The BLM scale (^) might not describe accu¬ 
rately the typical virtuality of the gluon in such a case. The possibility of a change 
of sign in p{s) is not excluded in [0. However, the examples given there suggest 
that the positivity of p{s) may indeed be generic. This is certainly the case for 
the vacuum-polarization D function: p{s) > 0 for any s, and for the Bjorken sum 
rule, as we show in section 4.2. In examples drawn from heavy-quark physics, p{s) 
becomes negative only as an artefact of the UV cut-off, and it may well be that p{s) 
would be non-negative for an appropriate choice of the UV regulator. 

In the rest of this paper we restrict our attention to cases where p{s) > 0 for any 
s. This will be a crucial assumption for most of our results, which can be checked 
only if p{s) is calculated exactly, or if suitable general theorems can be proven. 

We begin by showing that the resummed result in the large-/3o limit is a Ham¬ 
burger function of the renormalized coupling. Our starting point is expression ([TT| ) 
for the resummed result in the large /3o approximation. Since the integral is scale- 
and scheme-invariant, we can work in the V scheme where C = 0, and use the renor¬ 
malization scale = Q^, without loss of generality. In the large-/3o approximation, 
it is convenient to dehne 2 ; = as the coupling constant. Thus we get: 


/ OO 

P(s)- 

-OO 


^ -ds = r 


1 “h SZ J —oc 1 -|- SZ 

where (j){s) has the property that p{s) = d(l){s)/ds. Then we define the moments 


( 22 ) 


fi = (-1)*C,>'((?") = / s'd^(s) 


(23) 


for z > 0. 

If the moments are hnite and p{s) is a non-negative function, i.e., 0(s) is a non¬ 
decreasing function, then f{z) is a Hamburger function [^. Although f{z) has 


a cut on the real axis, it is well dehned for complex arguments. In the special case 
where the integration is limited to positive values of s, i.e., p{s) = 0 for s < 0, f{z) 
is a Stieltjes function, which is well dehned on the positive real axis, since it has 
a cut only for negative real .2 values. In this case, the moments fi are all positive 
and the perturbative series fi{—zy is a Borel-summable asymptotic series. 

However, QCD momentum distribution functions are not Stieltjes functions, since 
long-distance effects make p{s) non-zero for negative scales s. 

At the present stage of knowledge of QCD, the statement that Ares should be a 
Hamburger function for a generic physical observable is just an assumption, equiva¬ 
lent to assuming that p{s) is positive. If one knows only the hrst n perturbative coef- 
hcients, one cannot construct the exact distribution and therefore cannot prove that 
it is indeed a Hamburger function. However, Hamburger functions satisfy certain 
consistency conditions: these can be checked already using the hrst few calculated 
moments, and functions that do not satisfy them cannot he Hamburger functions. 
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One such criterion is the following: consider the determinants D{m,n) dehned by: 


D{m,n) = 


fm fm+1 

fm+1 fm+2 


fm+n fm+n+1 


fm+n 

fm+n+l 


fm+2n 


In the case of a Hamburger function, the determinants D{m, n) for even m and any 
n are positive [^. Another criterion is provided by the PA’s discussed shortly: if 
the poles of the z[N — 1/A^] PA’s include complex pairs, the momentum distribution 
function cannot be a Hamburger function. 


3.2 Properties of PA’s of Hamburger Functions 


We now use the characteristics of Hamburger functions to draw conclusions on the 
use of PA’s for resumming the perturbative series for a positive distribution func¬ 
tion. We hrst draw attention to the theorem that a PA of a Hamburger function 
has only real roots and positive weights. This means that successive PA’s dehne 
approximations to the distribution function in terms of 6 functions with the corre¬ 
sponding locations and weights. The hrst such approximation at NLO is provided 
by the x[0/l] PA, which is equivalent to the BLM method P], and higher-order 
approximations are provided by the x[N — 1/W] PA’s. These approximations to the 
momentum distribution function have the advantage of yielding scale- and scheme- 
invariant results when one integrates, in the large-/do approximation. As we discuss 
in more detail below, these approximations chosen by the PA are identical to the op¬ 
timal choices based on the all-order distribution function according to the Gaussian 
quadrature method for numerical integration. Then we show how, for a Hamburger 
function, the hrst few coefficients of the perturbative series can be used to set rig¬ 
orous bounds on the all-order momentum distribution function. Finally, we discuss 
the non-convergence of PA’s due to the presence of IR renormalons, and mention 
the non-physical case of a Borel-summable series that dehnes a Stieltjes function, 
for which PA’s do converge. 

In the language of momentum distribution functions, the perturbative coefficients 
may be interpreted as moments of the exact distribution function, as in (|^. For a 
Hamburger function, it can be shown that the z[N — 1/W] PA constructed from 
the partial sum zYl,i=o fii—zY, where n = 2N — 1, may be written as: 


N 


fiz) ~ z[N-l/N] = 


TiZ 


(24) 


where (i) the locations of the poles - 
Tj are positive for alH = 1, 2, ...W. || 


^^l + qiZ 

1/qi are on the real axis, and (ii) the weights 


II Note that the signs of the residues ri/qi are determined by the signs of the 
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3.3 Approximations to the Momentum Distribution Func¬ 
tion 


It is easy to construct approximations to p{s) which yield (pi] ) when substituted in 
(P^. These approximating distributions, which we denote by Pn{.s), are given by 


N 

Pn{s) = '^ri6{s - Qi). 
i=l 


(25) 


In words, the approximations are given by sums oi N 6 functions, located at the 
points Qi = —1/pi, where the z = pi are the PA pole locations, and with weights r,. 
Therefore PA’s have very natural descriptions as approximations to the momentum 
distribution function. Correspondingly, 0(s), the indefinite integral of p{s) is ap¬ 
proximated by monotonically non-decreasing piecewise-constant functions with N 
steps: 

N 

^n{s) = - Qi) (26) 

i=l 

where 6{s) is the Heaviside function. 

We note that the naive nth-order perturbation series has no natural description 
in terms of the distribution function. Formally, it may be expressed as a singular 
function composed of the n hrst derivatives of the Dirac 5 function: 

zY,Ckz'‘= Pis)——ds (27) 

fc=o i + sz 

for p{s) = Ylk=oi^k/k\)5^’^\s). It is clear that the corresponding 0(s) is in general 
not a good approximation to 4>{s): for one thing, it is not a monotonically non¬ 
decreasing function. 


3.4 Relation to the Gaussian Quadrature Method 

To underline further the utility of PA’s, we note that there is a formal mathematical 
relation between them and the Gaussian quadrature method for numerical inte¬ 
gration, which may provide an opportunity to extend the resummation to include 
non-leading terms in Pq. 

The basic quadrature problem is a generalization of ours, namely to hnd a for¬ 
mula for the numerical integration of a given arbitrary function y{s) with respect to 
a positive weight function p{s): 


rb N 

y{s)p{s)ds = y{si)wi + cn 

^ i=i 


(28) 


where cat is the error. One seeks the sampling points s* and weights Wi which 
minimize the error e^v. The well-known A^th-order Gaussian Quadrature method 


reviewed in [211 is to choose the N sampling points and corresponding weights such 
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that any polynomial function of maximal order n = 2N — 1 substituted for y{s) will 
be integrated exactly. The error for any other smooth function will in general be 
small, since it results only from the difference between the exact y{s) and its best 
weighted polynomial approximation. The condition that any polynomial of maximal 
order n will be integrated exactly is: 

/ s^p{s)ds = J2si'"u!i (29) 

2 “ i=i 


for k = 0,1, ...n. 

It is easy to see that the required sampling points and weights can be obtained 
from the a;[iV — 1/iV] PA of the following Hamburger function: 



N 


^ p{s)ds = t 


1 + zs 


t^^l + ZSi 


+ Oiz^^) 


(30) 


In other words, the Gaussian quadrature formula for numerical integration is ob¬ 
tained by replacing the exact continuous distribution function p(s) by a weighted 
sum of 6 functions, with locations and weights determined by the x[N — 1/N] PA 
of the corresponding Hamburger function. 

Specializing now to the physical QCD problem, we start with the following gen¬ 
eral expression for Ares- 

/ CO POO 

x{e'^^‘^Q‘^)p{s)ds = / x^{e^Q‘^)p{s)ds (31) 

-oo J—oo 

In order to relate our problem to the Gaussian quadrature integration problem, we 
identify the function y[s) over which one integrates in (^Sj) with the running cou¬ 
pling constant: y[s) = x ^and the weight function p{s) with the Hamburger 
momentum distribution function. The inhnite integration range is not expected to 
cause any trouble, since p{s), being renormalized and IR hnite, is expected to vanish 
for large positive and negative arguments. Thus our physical resummation problem 
is very close to the Gaussian quadrature integration problem described above. The 
one difference in our physical problem is that the running coupling constant is in¬ 
tegrated with respect to a weight function that is not fully known. The only pieces 
of information we have about this function are its hrst few moments Cq through 
Cn- This is to be contrasted with the mathematical integration problem, where 
the weight function p{s) is known exactly, but we limit the numerical calculation 
to a certain order, for other reasons. However, at any given order, the choices of 
scales and weights furnished by the x[N — 1/N] PA of the corresponding Hamburger 
function is identical with the optimal choice, according to the Gaussian quadrature 
formula based on the exact weight function. Specihcally, we obtain 


Ares ^^TiX^ 
i=l 


(32) 
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where the locations g* and weights r* are computed from the z[N — 1/iV] PA (|2J) 
of the leading order of the large-/3o series. As already remarked, the usual BLM 
approach corresponds to A^ = 1. Since the running coupling cannot be described as 
a hnite-order polynomial in the scale, but rather as an inhnite-order series, the PA’s 
cannot yield the full all-order resummation, but they do yield ‘optimal’ approxima¬ 
tions to it. 


One can also use equation (|3^) to resum effects that are non-leading in /do in 
the running of the coupling constant, simply by using a higher-order formula for 
involving /3i, /32 and so on. We note that this procedure for using PA’s to 
resum the series is different from simply constructing the PA’s of the partial sums 
that contain non-leading terms in It seems likely to give more precise results, 
though this requires further study. 


3.5 Bounding the Momentum Distribution Function 

In practice, one wishes to use the PA method and the momentum distribution for¬ 
malism for QCD observables for which we know only a few low-order perturbative 
coefficients. As we have shown, PA’s provide an approximation for the distribution 
function, and one then integrates over the coupling constant with this approximate 
distribution function as a weight. Mathematically, this is analogous to a moment 
problem, namely the construction of a distribution function from its moments. If 
we assume that the momentum distribution is non-negative, we obtain the so-called 
Hamburger moment problem. We have already seen how the x[N — 1/N] PA’s are 
related to approximations to the integral of the momentum distribution function 
involving N steps In this subsection, we use some further mathematical the¬ 

orems related to the Hamburger moment problem to show how the all-order 
distribution function can be bounded rigorously - assuming that it is non-negative - 
by continuous upper and lower bounds that are based on PA’s, using only the hrst 
few coefficients of the perturbative series. 

We start with the dehnition (^) of the Hamburger function, and assume that 
the moments fi dehned by (|^) are known for i = 0,1...2M (M > 1). As before, we 
consider the partial sum zJffl^Qfi{—zy: n = 2M + 1, and construct the diagonal 
PA as usual: 

A^{z) 


z[M/M + 1] = z- 


(33) 


where A^{z) and B^~^^{z) are polynomials of orders M and M -|- 1, respectively. 
In principle, we cannot construct this PA if / 2 M +1 is unknown. However, the value 
of / 2 M +1 will eventually be eliminated from our hnal results for the bounds on the 
distribution function, so knowing it is actually not essential for our present purpose. 

Following the method described in part H, section 3.2 of [^, we construct an¬ 
other PA: 

fiMf \ 

z[M/M] = z^^^^ (34) 


D^{z) 
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where C^{z) and D^{z) are both polynomials of order M. We note that z[M/M] is 
an off-diagonal approximant, and is therefore not invariant nnder scale and scheme 
transformations, in contrast to the z[M/M -|- 1] PA An important observation is 
that z[M/M], like the z[M/M -|- 1] PA, has only real roots. 

Using the two PA’s (|33| , we introdnce an anxiliary variable w and dehne 

g{z) = D^{z)B^+^{w) - D^{w)B^+^{z) (35) 


and 


h{z) = C^{z)B^+^{w) - D^{w)A^{z) 


(36) 


It is helpfnl to think of g{z) and h{z) as polynomials in z, with tc-dependent co¬ 
efficients. One can then develop some intnition for several interesting resnlts on 
the ratio h{z)/g{z), provided in ref. |^. We do not prove them here: rather, we 


snmmarize 


them briefly and refer the interested reader to 0 for details: 


The qnantity zh{z)/g{z) is an approximation to the original Hambnrger fnnc- 
tion f{z), with the property that: 


(J(z)/z)a(z)-h(z)=0{z^^+') 


(37) 


• The qnantity zh{z)/g{z) can be written in the form 

h{z) zpj 


(38) 


where pj and Q depend on w. In particnlar, = —1/w for some k. 


Mimicking the relationship (^) between f{z) and 0(s), one can constrnct from 
(^) a monotonically increasing integral distribntion fnnction 'ijj{s), 




M+l 

- 


i=i 


0 ) 


(39) 


that, becanse of the property mentioned above, has a point of increase at 
s = —1/w. 


• One can then dehne 


V'-(s) = lim V'(s) (40) 

w=(—l/s) 

Ms)= lim U(s) 

W = [ — l/s) + 

and it can be shown that, for any s, 


V'_(s) < 0(s) < V’+(s) (41) 


**However, our final bounds are scheme and scale independent. 
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These results imply that the assumption that the distribution function p{s) is non¬ 
negative allows one not only to extract some good approximations to the all-order 
distribution from the hrst few coefficients, but also to evaluate the errors, by the 
construction of rigorous bounds on the distribution. The bounds (^T]) , just like the 
x[N — 1/N] PA’s are renormalization-scale and scheme invariant. This is a direct 
consequence of the fact that zh{z)/g{z) is a z[M/M -|- 1] PA for any value of the 
parameter w. 

As an illustration, we present here the explicit results for the bounding distri¬ 
bution function in the case where only three coefficient in the perturbative series 
are known, namely M = 1 in the above general analysis. This is both the simplest 
case where such bounds can be formed, and the maximal order for which QCD 
observables have been fully calculated up to now. Thus, we start with a series: 


f{x) = x{Cq + CiX -1- C 2 X^) 

and use the a;[1/1] and a:[1/2] PA’s to obtain the following bounds: 


V’-(s) = 


((Tog -|- CiY 
Cos2 + 2Cis + C 2 


0 s + 


Cis -\- C2 

CoS + Cl 


and 


Ms) = + 


C'oCa - Cf 


(42) 


(43) 


(44) 


+ 2Cis + C 2 

We see that the bounds approach each other as the series resembles more closely a 
geometrical series, as intuitively expected. 

Finally, we note a certain complementarity between this method of constructing 
bounds for the distribution function, and the previous PA estimates of the distribu¬ 
tion function. This arises because the bounds are based on the perturbative series 
zJ27=ofii~^)^ evaluated to some even order n = 2M, whereas the a;[A^ — 1/A^] PA 
requires the knowledge of the perturbative coefficients to some odd order, n = 2A^—1. 


3.6 Convergence of Higher-Order PA’s 

Knowing that PA’s resum a generic perturbative QCD series, one may naively expect 
that they should converge at high orders. However, it is well known that perturba¬ 
tive series in QCD do not contain full information about long-distance effects. In 
particular, IR renormalons appear in the formalism of the momentum distribution 
function through the combination of the non-vanishing of the distribution function 
p{s) at negative s with the Landau pole in the running of the coupling constant seen 
in Mathematically, this implies that A^es has a cut on the real axis. From the 
physical point of view, this can only be avoided |^, by postulating freezing of 
the coupling constant 


23|: for a recent discussion, see p4 


From a mathematical point of view, it is interesting to look at another limit in 
which Aj-es is well dehned on the positive real axis, namely when p(s) vanishes for 
negative s. If one assumes that p(s) = 0 for negative s in (^), in addition to the 
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assumption that p{s) is non-negative, as already mentioned one obtains a Stieltjes 
function: _ 

rather than a Hamburger function. In this case, f{z) is well dehned on the positive 
real axis, though it still has a cut on the negative real axis, and therefore the formal 
Taylor expansion of f{z) ~ fi{—zy has a zero radius of convergence. Unlike 

the Hamburger series, where the coefficients have no dehnite sign, here the series 
oscillates in sign, since the moments ft are all positive. There is a theorem for 
Stieltjes series that higher-order PA’s converge to the true Borel sum of the series, 
even though the perturbative power series diverges. 

In QCD, both the UV and the IR parts of the momentum distribution function 
are expected to exist, so we do not have the case of a Stieltjes series. Moreover, 
studies of different examples indicate that the errors of PA’s are particularly large 
when the IR and UV parts are about equally important. However, we should like 
to stress that the non-convergence of increasing-order PA’s is expected to be much 
softer than that of the corresponding partial sums in perturbative QCD. Whilst PA’s 
oscillate around the Cauchy Principal Value of the Borel-resummation integral |^, 
the partial sums diverge badly due to the zero radius of convergence of the series. 
Thus the effects as well as the reasons for the divergences are different. 


4 Some Worked Examples 

4.1 The Vacuum Polarization D Function 


As a hrst example of the general results of Sec. 3, we choose the particular case of 
the vacuum-polarization D function 


D{Q^) = 


dn(g2) 

dQ^ 




/ 


1 + A£)(g^) 


(46) 


to test our method. We neglect here a small light-by-light contribution. For this 
example, the all-order resummed result is known in the large-/?o approximation [p!8| . 

Following 


13, 20 : see also 25 


we write the momentum distribution function 


in the form 


Pd{s) 


= 8C/e" X 

I e® -1- (1 -F e^) [L 2 (-e®) + sln(l -h e^)]| 


1 + S + — 



-F(l-FeQ L 2 {-e Q-sln(l-Fe Q 


(47) 
s < 0 

s > 0 


where L 2 {x) = — y ln(l — y) and Pd{s) is the weight function for the 
tion integral: 



Pd{s) — 


x^{Q^) 

(s + C)/3ox«(g2) 


resumma- 

(48) 


16 























where we choose the renormalization scale jj? = but still allow for an arbitrary 
scheme R. Using equation (|i8|) we can obtain any required coefficient d^{Q‘^): 

/ OO 

Pd{s) {s + Cy ds (49) 

-OO 

and the nth-order partial sum is: 

n 

A„{Q^) = 4{fkx’‘(Q^)Y (50) 

i=0 

whose accuracy we now compare with PA’s. 

Fig. m presents the partial sums (0), as a function of n for = 2 GeV^. 

At this low value of the perturbative series starts to diverge already at relatively 
low order: n ~ 5, so the differences between various renormalization schemes and 
other calculational approaches is readily apparent. The horizontal continuous line 
in Fig. 1^ represents the all-order resummation of the leading Pq terms, as calculated 
by taking the Cauchy Principal Value of the integral in (^8|). The two horizontal 
dash-dotted lines correspond to the maximal uncertainty which is inherent to the 
integration in (^81) , due to the hrst IR renormalon. 

Superposed on the all-order resummation results in Fig. we show the naive 
perturbation theory partial sums in various schemes with . In this case, the 

MS partial sums of increasing order converge quite nicely to the all-order result, 
whilst the the U-scheme results are totally misleading. It is important, however, 
to realize this relative success of MS with p? = in this case has no known 
theoretical basis. There are other examples where MS with is not a 

good choice, such as the Bjorken sum rule considered in [^. It is clear that, if one 
chooses to evaluate partial sums, it is essential to have some criterion for choosing 
an appropriate renormalization scale and scheme |^. 

An example of a judicious choice of scale is provided by the BLM criterion, which 
sets the scale such that the NLO contribution vanishes, as seen in Fig. |I], where the 
leading-order and NLO BLM results are indeed the same. The alternative and 
generalization to higher orders that we suggest, namely the x[N — 1/V] PA, does 
not have any scale ambiguity. The x[N — 1/N] results presented in the hgure are on 
par with what one gets after optimal tuning of renormalization parameters within 
the usual schemes, and certainly much better than the results with a generic choice 
of renormalization scheme and scale. 

Next we consider the equivalence between taking the x[V — 1/V] PA of the 
perturbative series and approximating the momentum distribution function by a 
weighted sum of 6 functions, as discussed in section 3.3. The large-/3o D-function 
momentum distribution of (|48|) is indeed non-negative, and the resummed result 
Ad{Q^) is therefore a Hamburger function. Hence the results of section 3 are fully 

As mentioned in section 2, the scale and scheme parameters play the same mathematical role 
in the large-/3o limit, so there is only one parameter to tune in this case. Beyond the large-/3o 
approximation, there are more parameters to specify. 
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applicable to this physical example. Since the hnal result in our method is scheme- 
and scale-independent, we are free to choose the Id-scheme with = Q^, which 
simplihes the calculations. Starting with the nth-order perturbative series, where 


n = 2N -1: 


Ad ^ {Q^)Y,dY{Pox^{Q^)y = 


i=0 


fdo i=o 


(51) 


with the coefficients df calculated from with C = 0, we construct the 2 ;[A^ —1/iV] 
PA. We determine numerically the locations of its poles and their corresponding 


residues. As guaranteed by the general theorem |^|, all the PA poles —1/di are real 
and all the weights r* are positive for any i. We then hnd, for every z[N — 1/A^] PA, 
the corresponding approximations to the momentum distribution function Pd{s)'- 


N 


Pn^s) = ^ri(5(s-di). 
and to its indehnite integral (f>D{s) 

-?j) 


(52) 


2=1 


N 


(53) 


2 = 1 


as described in section 3. 

We illustrate in Fig. the way in which the momentum distribution function p{s) 
for the D function is approximated by a sum of A^ <5 functions, corresponding in panel 
(a) to the x[N — 1/N] PA’s for A^ = 1 through 3, and in panel (b) to the higher- 
order case N = 12. In each panel, the exact continuous distribution Pd(s) of (|i5|) 
is shown here as a solid line. Superposed on it, we plot the locations and weights 
Tj of the 6 functions which compose the function p^{s) of (^). At leading order, 
the x[0/l] PA coincides with the BLM method, as can be verihed from the value of 
the scale parameter: ~ 0.975, which corresponds to — 1-628a/^. The 

convergence of higher-order PA’s is particularly clear in panel (b) of Fig. Q - note 
that the vertical scale is logarithmic. 

Fig. 0 shows the corresponding integral 0(s) of the distribution function for the 
D function, including both the exact distribution (continuous line) and the 

hrst three approximations of (^Bj), that correspond to the x[N — 1/N] PA’s for 
A^ = 1, 2, 3. We see how the N steps imitate the shape of the all-order momentum 
distribution, based only on knowledge of the hrst few moments of the distribution. 

Finally, in Fig. || we illustrate the PA-based bounding technique of section 3.5, 
for the same example of the vacuum-polarization D function. The exact integral 
distribution function </>d(s) is again plotted as a continuous line. We have calculated 
the upper and lower bounds in the two simplest cases, the hrst being based on the 
a;[1/2] and a;[l/l] PA’s and requiring knowledge of do, di and d 2 (plotted as a dashed 
line), and the second being based on the a;[2/3] and a;[2/2] PA’s and requiring the 
knowledge of do through d^ (dotted line). As before, the calculation was done in the 
V scheme. Just as for the x[N — 1 /N] PA, if one uses the bounds on the distribution 
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function in order to integrate over the coupling constant, one obtains a scale- and 
scheme-independent result. 

We conclude this subsection with an interesting empirical hnding, for which we 
have no good explanation. As can be seen in Fig. ^ the arithmetic mean of the 
upper and lower bounds on the distribution function is quite close to the exact all¬ 
order result. This may very well be a coincidence, but also might have some deeper 
theoretical justihcation. 


4.2 The Bjorken Sum Rule 


As an indication that the above example is not isolated, we now consider the per¬ 
turbative QCD series for the Bjorken sum rule, again in the large-/3o approximation. 
This series is known to have a particularly simple structure in the Borel plane 1^^ : 


B{u) = 


8 


9(1 + m) 18(l + iM) 9(1 -m) 18(l-iM) 


(54) 


where B{u) is the Borel transform, containing only four simple poles. Using an 
inverse Laplace integral [jT^ we can derive from (|5^) the corresponding momentum 
distribution function: 


p{s) = 4 


,-2s 


0(S) + 


/16 

\ 9 ^ 


10 


2s 


e(-s) 


(55) 


This momentum distribution function is plotted in Fig. 5, where we see explic¬ 
itly that it is positive, and hence dehnes a Hamburger function, as is the vacuum- 
polarization D function in the large-/3o limit. 

It is clear that the steps carried out for the previous example, namely the approx¬ 
imation of p{s) (|55[) by a sum of 6 functions, the evaluation of the corresponding 
integral 0(s) and the establishment of PA-based bounds, can be carried out in a 
similar way, but we do not enter here into the details. 


5 Conclusions 

This paper has been devoted to analyzing the PA method in QCD, and to under¬ 
standing the reasons for its success in resumming perturbative series, which were 
previously unclear B i- We now understand that rigorous conclusions can be 
drawn in the large-/5o limit regarding the use of PA’s. In particular, we hnd that 
the x[N — 1/N] PA’s are the most appropriate for resumming the series, because of 
two important characteristics. 

• They are scale and scheme invariant [^. 

• If the momentum distribution of a virtual gluon in a generic QCD observable 
is non-negative, the resummation integral - which is a weighted average of the 
running coupling - dehnes a Hamburger function. In this case, the x[A^ — 1/A^] 
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PA is also the result of an exact integration of the coupling constant over 
the A^th-order optimal approximation to the momentum distribution function. 
Therefore, the resulting resummation makes full use of the hrst n + 1 = 2N 
coefficients of the series that are known. 


We also saw how the assumption of positivity of the momentum distribution func¬ 
tion allows one to construct scale- and scheme-invariant bounds on the all-order 
distribution, using knowledge of the hrst few perturbative coefficients. 

There are two main questions that our work raises. The hrst is: What are the 
conditions for our conjecture on the positivity of the distribution function to be 
valid? In some cases, where the large-/3o all-order resummation has been performed, 
this can be checked explicitly, as was done for for the D function in [|1^ and for the 


Bjorken sum rule in this paper. In other cases, where only a few hrst coefficients are 
known, it would only be possible to disprove the conjecture, as discussed in section 
3.1. It would be very interesting to hnd theoretical justihcation why the momentum 
distribution function should be positive. 

The second question is: How can one extend the PA’s method, and use it outside 
the large-/do limit? It is obvious that the rigorous results we have obtained in this 
limit cannot be extended in a straightforward manner. The naive approach of using 
the x[N — 1/N] PA of the full series is not very well motivated by the momentum 
distribution concept, and gives a wrong functional dependence on the number of 
colors and light havors. An alternative, which is motivated by the Gaussian quadra¬ 
ture integration procedure, is described in section 3.2. The basic idea is to use the 
scales and weights of the corresponding large-/3o series, and a higher-order formula 
for the running coupling constant. This method does not fully use the perturbative 
coefficients that are known, so there is still room for further improvement. 

Despite the persistence of these open questions, we feel that the analysis of this 
paper has contributed to a useful theoretical foundation for the use of PA’s in ap¬ 
plications to perturbative QCD, and also clarihed in a useful way the relation of 
the BLM method to the PA approach. We hope that this paper may serve as a 
helpful building block in the search for an eventual optimal strategy for exploiting 
the information contained in perturbative QCD series. Higher-order QCD calcu¬ 
lations have recently taken several impressive steps forward. However, it seems in 
many cases unlikely that the following terms in the series will be forthcoming in 
the foreseeable future. Moreover, we shall in any case only have access to a hnite 
set of exact terms in any perturbative QCD series. Therefore tools to optimize the 
return on the investment made in exact calculations will always be a welcome input 
to making precision tests of QCD and understanding its place in an eventual unihed 
theory. We believe that both PA’s and the BLM method both contribute to the 
fashioning of these tools. 
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Figure 1: Increasing-order results for the vacuum-polarization D function in the 
large-/3o limit, as calculated in different renormalization schemes/scales - MS with 
jj? = Q^, the V scheme with = Q^, the BLM method, and diagonal PA’s. The 
horizontal continuous line corresponds to the Cauchy Principal value of the Borel 
integral of the perturbative series, and the horizontal dash-dotted lines describe the 
maximal intrinsic ambiguity of this integration due to the hrst IR renormalon pole. 
Note that the diagonal PA’s can be constructed starting with any dehnition of the 
coupling, i.e., using any scale and scheme, to give the results that are plotted. Note 
also that leading-order result in the BLM method coincides exactly with the a;[0/1] 
PA. 
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Figure 2: The momentum distribution function p{s) in the large-/3o limit, for a 
virtual gluon in the vacuum-polarization D function. In each panel, the solid line is 
the exact result for p{s), and the different symbols correspond to the locations and 
relative strengths (weights) of Dirac 6 functions, as determined by diagonal PA’s. 
As discussed in the text, the x[N — 1/N] PA chooses the locations and weights 
such that the 2N hrst moments of any function integrated with respect to p{s) are 
reproduced exactly. Panel (a) shows low-order PA’s on a linear scale, and panel 
(b) shows a representative high-order PA on a logarithmic vertical scale and an 
expanded horizontal scale. The convergence of the PA’s to the true momentum 
distribution function is clearly visible. 
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Figure 3: The integral of the momentum distribution function in the large (3q limit, 
for a virtual gluon in the vacuum-polarization D function. The continuous line 
represents the exact result for (j)£,{s), and the other lines describe different approxi¬ 
mations to 4>d{s) that correspond to the x[N — 1/iV] PA’s of the perturbative series 
for N = 1,2,3. The A^th approximation to 4>d{s) is a piecewise-constant function 
composed of N steps, with heights determined by the weights of the 6 functions 
provided by the corresponding PA’s. 
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Figure 4: The integral momentum distribution function in the large (3q limit, for a 
virtual gluon in the vacuum-polarization D function. The continuous line represents 
the exact result for 0 _d('S), and the other lines describe the upper and lower bounds, 
as well as their averages, that one can construct as described in section 3.5 from 
the first three (dashed lines) and hve (dotted lines) coefficients of the perturbative 
series. Note that the average is very close to the exact distribution function. 
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